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Abstract. We study the domain of validity of Perturba- 
tion Theory (PT), by comparing its predictions for the 
reduced skewness, S3, and kurtosis, S4, of the projected 
cosmological density field, with the results of N-body sim- 
ulations. We investigate models with different linear power 
spectra and consider as physical applications both angu- 
lar galaxy catalogues and weak lensing surveys. We first 
find that the small-angle approximation for the predicted 
skewness provides a good match to the exact numerical 
PT results. On the other hand, results from non-linear 
simulated catalogues agree well with PT results on quasi- 
linear angular scales, which correspond to scales larger 
than about 1 deg in the applications we have considered. 
We also point out that, on smaller scale, the projection 
effects tend to attenuate the effects of the strong nonlin- 
earities in the angular skewness and kurtosis. 

Key words: Cosmology: Dark Matter, Large-Scale Struc- 
tures, Gravitational Lensing 



1. Introduction 

One of the principal goal of the large galaxy surveys that 
should be available in the near future (SDSS, 2DF, deep 
surveys with the MEGACAM project, ...) is the accurate 
determination of the power spectrum of density fluctu- 
ations, P(k). Its interpretation in terms of cosmological 
models, however, requires an understanding of the way 
galaxies trace the underlying matter distribution. This is 
a general problem that is bound to become even more 
crucial with these new observational data. 

It has been stressed recently (Fry & Gaztanaga 1993, 
Gaztanaga & Frieman 1994, Fry 1996, Bernardeau 1995, 
hereafter B95) that one way to address this problem is 



to consider higher-order correlation functions or higher- 
order moments of the local density probability distribution 
function. This approach is particularly attractive because, 
in the case of Gaussian initial conditions, Perturbation 
Theory (PT) provides precise quantitative predictions. It 
has now been well established that the p-order cumulants 
of the local density field ($ p ) c are expected to behave as 



2\P-2 



(1) 
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on large scales (Fry 1984, Goroff et al. 1984, Bouchet et 
al. 1992, Bernardeau 1992). The S p parameters, which 
quantify the departure from Gaussian behavior, depend 
however on the window function applied to the field. In 
Bernardeau (1994) a prescription is given for the PT cal- 
culation of all these coefficients for a 3D top hat window 
function. These results have been subsequently checked in 
detail and found to be very accurate when compared to 
numerical simulations (Baugh, Gaztanaga & Efstathiou 
1994, Baugh & Gaztanaga 1995). Perturbation results for 
the 3D Gaussian window function (Goroff et al. 1986, 
Juszkiewicz, Bouchet & Colombi 1993) for S3 and S4 
have also been successfully tested against numerical sim- 
ulations. Another fruitful direction of investigation is the 
calculation of the high-order correlation functions for the 
projected density (B95). This is an interesting domain of 
investigation since angular galaxy catalogs contain more 
objects and volume than 3D catalogs and therefore allow 
the determination of many more parameters of the hier- 
archy (1). 

There is also a new mean of investigation, which is 
still in an embryonic state but might reveal extremely 
fruitful, for the projected density. The measurement of 
the gravitational weak shear induced by the large scale 
structures in deep galaxy catalogs allow in principle to 
have access to the correlation properties of the projected 
mass. The resulting polarization maps could allow the de- 
termination of these correlation functions at the level of 
the two-point function (Blandford et al. 1991, Miralda- 
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Escude, 1991, Villumsen 1996, Jain & Seljak 1997, Kaiser 
1996) or even for higher orders (Bernardeau, van Wear- 
beke & Mellier 1997). 

Although the projected density is obtained with com- 
pletely different methods in galaxy catalogs or in weak 
lensing surveys, the properties of the projected density 
as yielded by Perturbation Theory are addressed in very 
similar ways. These two cases differ only in the shape 
of the required selection function, whereas identical ap- 
proximations are made in the course of the calculations 
for the derivation of the analytic expressions. We recall 
here that Perturbation Theory results are valid at large 
scales, where the variance is small. For 3D filtering, in- 
vestigations in numerical simulations have proved, for 3D 
filtering, that PT results are valid for scales above about 
10 ft- _1 Mpc, where the variance of fluctuations approaches 
unity. In the case of the projected density, however, a given 
angular scale cannot be straightforwardly associated with 
a physical scale, since it corresponds to a superposition 
of different scales. A priori it is therefore difficult to as- 
sess the validity domain of PT results, even in light of 
the 3D cases. Another concern is the use of the small an- 
gle approximation. This is a mathematical approximation 
that allows one to dramatically simplify the calculations. 
So far it has not been possible to get closed analytic for- 
mulae without its use. This approximation might be not 
fully valid when the smoothing angle is above 1 degree 
(see B95). 

The aim of this paper is therefore to investigate the 
validity of both the small angle and the PT approxima- 
tions. Numerical simulations with iV-body codes are used 
to check the validity of PT at small scales, and to give 
better clues to the importance of the fully non-linear cor- 
rections that may affect the quantities we are interested 
in. The effect of the small- angle approximation is more 
particularly investigated by direct Monte-Carlo integra- 
tion of the S3 coefficient. In this case, the problem is not 
the validity of the Perturbation Theory approach, but the 
accuracy of the mathematical approximations that were 
made in the course of the calculation. 

In §2, we specify the models that are used to illustrate 
those calculations. We also give the expected values for S3 
in those cases. In §3 and §4, we confront those predictions 
with the Monte-Carlo calculations and the simulations. 



2. The Physical Models 

To illustrate our calculations we consider two different 
cases for the linear density power spectrum. We also 
consider two different physical applications: the selection 
function corresponding to angular catalogs and that cor- 
responding to detection of weak shear. 



2.1. The power spectra 

In our first investigation we use the power spectrum de- 
rived from the angular APM Galaxy Survey (Maddox et 
al. 1990), which has been found to be well described by 
(Baugh & Gaztanaga, 1996) 



PAPM(k) oc 



with 



[1 + (fc/fcc) 
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In the second case, we use the standard CDM model 
with fibaryon = 0.05, H = 50 fcm/s/Mpc n = 1, A = 0, 
and an initial Harrison-Zcl'dovich spectrum. For conve- 
nience the mass fluctuation power spectrum is approxi- 
mated by a simple analytic fit (similar to the ones pro- 
posed by Bond & Efstathiou 1984), 



PcDAt(k) = A- 
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with u = 1.13, a = 



6.5 



1.7 



3000 r " 3000 r ~ 3000 r r 

0.5, where the dimensional quantities have been expressed 
in such a way that c = Hq = 1. 

2.2. The angular galaxy catalogue selection junction 

In this case the local projected density field ui in the di- 
rection 7 is given by 



dr r 2 F(r) S(r, 7), 



(5) 



where 5(r, 7) is the local over-density in the direction 7 
at the distance r, and F(r) is the normalized selection 
function for the catalogue, 



dr r 2 F(r) = 1. 



(6) 



In case of the APM angular survey the radial selection 
function is given in Gaztanaga & Baugh (1997, hereafter 
GB97). This function can be approximated by: 



F{r) cx r~ h exp[-r^/£> 2 ]. 



(7) 



We will use b ~ 0.1 and D ~ 335ft. -1 Mpc, which provides 
a good fit to the APM selection function, as shown in 
Figure H (see §|], below). 

We can calculate the moments of the distribution of 
the local projected density Sg, smoothed at scale 8 (see 
Bernardeau 1995, Polio & Juszkiewicz 1997) using Per- 
turbation Theory and the small angle approximation. The 
second moment, or variance of angular counts-in-cells, at 
the smoothing scale 6, defined in (20), is given by 

(S 2 e ) = [dr r 4 F 2 {r) f kdk P(k) W 2D {k9) (8) 
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Fig. 1. Comparison the theoretical (lines) and measured 
counts (histogram) in radial shells for a mock catalogues. The 
dashed line is a model while the continuous line corresponds 
to the input APM selection function. 



with 



R* = 



J r 8-2(n+3) dr pi3( r 



(12) 



rj r 5-(n+3) d r F 2 (r)\ 

for a normalized selection function. The coefficient R3 is 
found in practice to be of order unity and to be very 
weakly dependent on the adopted shape for the selection 
function. Note however that the redshift evolution of the 
fluctuations has not been taken into account in this re- 
lation. This evolution should be taken into account for 
catalogs having a large depth. In this case, the geometri- 
cal factors for non-flat universes are also important (see 
G94). 

For the selection function given in equation (pfi we can 
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calculate s2 aL explicitly, 
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For b — and n = 0we find — —2= ~ 1.54, while for 

13 3^3 ' 

6 = and n = — 1, closer to the APM case, R$ 
1.21, comparable to the values given in G94. 
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where W is the window function in k space. This variance 
can also be related to an area average over the two-point 
angular correlation (e.g. Gaztanaga 1994, hereafter G94). 
Note that this expression does not rely on perturbative 
calculation if the power spectrum that is used is the final 
one. It uses, however the small angle approximation. In 
the following we will use exclusively the top-hat window 
function so that 
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where J\ is the spherical Bessel Function. The third mo- 
ment of smoothed angular fluctuations, defined in (21), is 
given by 



<^ 3 > = 
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2.3. The weak lensing efficiency function 

In the case of weak lensing the local projected density can- 
not be directly observed. We can however observe the local 
convergence (see Kaiser 1995, Bernardeau et al. 1997) k, 



kdkW 2 2 D (k9)P(k)x 
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So that, in case of a power law spectrum P(k) ~ k n , 
we have (B95), 
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where the efficiency function, E(r), depends on the red- 
shift distribution of the lenses. For simplicity we assume 
that the lenses are all at redshift z = 1 and that we live 
in an Einstein-de Sitter Universe. In this case we have 

*> = ^(^) 2 . 

with 

r s = 2 - V2. 

Here the second and third moment are given by 



(16) 
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and 
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These moments can be calculated explicitly only for 
particular power spectra. In the case of a power law spec- 
trum of index n, we have (see Bernardeau et al. 1997), 
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In the following we will compare these P.T. results with 
results of numerical simulations. 

3. A Monte- Carlo integration for s 3 

The aim of these Monte-Carlo integrations is to check the 
validity of the small angle approximation in PT. The in- 
tegration is made here in real space by throwing points at 
random in conical cells and by calculating the averages of 
the second and third correlation functions. We then iden- 
tify the moments with those geometrical averages, 
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where £ 2 and £3 are respectively the two and three-point 
correlation functions in real space and V conc is the 'volume' 
of the cone of angle 8 and a radial distribution given by the 
selection function F(r). In Perturbation Theory they can 
both be expressed in term of the linear power spectrum, 



6(ri,r 2 ) = 



and 
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As shown in the appendix it is actually more convenient to 
express these quantities through a set of real space func- 
tions. Indeed, we can define the function ip(r) by, 



ip(r) 



k 2 dk P{k) sin(fcr) 



2tt 2 



(24) 



from the derivative of which the three point function can 
be expressed. The Monte-Carlo integrations can then be 
done in real space. They reduce to 6 dimensional integrals 




Fig. 2. Angular skewness S3 for galaxy catalogue selection 
function from Monte-Carlo integration (symbols) compared 
with small angle approximation PT predictions (lines). The 
open figures and dashed line include (linear) redshift evolu- 
tion. 



(see appendix) which can be done on common work sta- 
tions. 

In Figures || and ||, we compare the results of the 
Monte-Carlo integration (symbols) for S3 with the small 
angle approximation (curves) in PT. We show both the 
case with (linear) redshift evolution (open figures and 
dashed line), which would correspond to real observations, 
and also the case without redshift evolution, as corre- 
sponds to the comoving output from N-body simulations 
(see below). 

There is an excellent agreement in the comparison spe- 
cially in the weak lensing case. For the galaxy selection 
function there are some small discrepancies at the largest 
scales, but they are smaller than the typical errors in the 
estimation from the observations. Thus the small angle ap- 
proximation seems to be quite good for most applications 
up to 10 degree scales. 



4. The numerical iV-body results 

We study two shapes of the power spectrum: the APM 
model and the CDM model, as described in Section §2. 
We use several sets of N-body simulations with parame- 
ters shown in Table [l]. The letters refer to different real- 
izations of each model: 5 for APM1, 2 for APM2 and 5 
for CDM. The APM like models are described in detailed 
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Fig. 3. Angular skewness S3 for weak lensing efficiency func- 
tion from Monte-Carlo integration (symbols) compared with 
small angle approximation PT predictions (lines). The opened 
figures and dashed line include (linear) redshift evolution. 



in GB97, while the CDM simulations are from GB95, and 
correspond to the standard model: T = 0.5. 



Table 1. Simulation parameters 



run 


number 
of particles 


mesh 


(fc^Mpc) 


APMl(a)-(e) 


126 3 


128 3 


400 


APM2(a)-(b) 


200 3 


128 3 


600 


CDM(a)-(e) 


126 3 


128 3 


378 



For the APM like models we use the output time cor- 
responding to the measured APM amplitude, <j% — 0.8, 
before it is scaled up to account for clustering evolution 
(see Gaztahaga 1995). For CDM models we use erg = 1, 
unless otherwise stated. As we are using outputs in co- 
moving coordinates, there is no redshift evolution within 
one output. This is not the case in the real Universe, but 
we take this into accunt in the theoretical predictions in 
a straightforward way (G94, see Figures || and ||). 

To produce the angular catalogues we first select an 
arbitrary point in the simulated box to be the local 'ob- 



server'. We include a simulated particle at comoving coor- 
dinate r from the observer with probability given by the 
selection function F(r), in the case of galaxy clustering, 
or the efficiency function in the case of weak lensing, as 
described in Section §2. As the simulation is done in a peri- 
odic box, we replicate the box to cover the total radial ex- 
tent of the APM (over 1800/i _1 Mpc). The main difference 
with GB97, is that they apply the APM angular survey 
mask, including plate shapes and holes, whereas here we 
use full sky maps, which cover a larger area and have no 
boundaries. In our full sky mock catalogues there could be 
a certain degree of repetition, but the errors are estimated 
form the dispersion in different catalogues. By comparing 
the results from different box sizes we have verified that 
this replication of the box does not introduce any spurious 
effects. The total number of particles in the mock angular 
catalogues is about 9 x 10 6 for the APM galaxy selection 
function and about 25 x 10 6 for the weak lensing case. 

Figure |l] shows a comparison between the expected 
number of galaxies, N(r)r 2 Ar, at different radial depths 
(in comoving coordinates r) given by the input selection 
function, as compared to the measured counts in a mock 
catalogue. The dashed line is the model in equation (0) for 
b ~ 0.1 and D ~ 335/i _1 Mpc, while the continuous line is 
the input APM selection function (GB97). 

From each realization we produce several mock cat- 
alogues by choosing different positions for the observer. 
Because of the selection function, catalogues from differ- 
ent observers are not necessarily correlated. Three differ- 
ent observers are used for the smaller boxes. For the larger 
boxes the corresponding number is 10 observers. For these 
numbers, results from different observers do not seem to 
be correlated. Thus the total number of mock catalogues 
is 15 for APM1 and CDM, and 20 for APM2. 

Figures |I] shows the results for S3 and S4 for the galaxy 
CDM and APM2 mock catalogues. Results for APM1 cat- 
alogues, which are not shown in the Figures, agree well 
with the ones in APM2, indicating that the size of the 
simulation box is large enough. 

There is very good agreement between the N-body re- 
sults and PT theory at scales 6 > 2 deg for the APM 
model. For CDM, there is agreement within the errors for 
9 > 1 deg, although this is not as good as for the APM 
model. This could be due to finite volume effects, which 
seem more important for CDM (see below); note also that 
the CDM errors are larger and that the CDM maps come 
from a simulation with a smaller volume. Scales 6 ~ 1 deg 
in the CDM models could also be affected by projection 
effects that tend to overcompensate non-linearities (see be- 
low); this is more important for CDM which has a higher 
normalization (erg = 1). 

Figure || shows a comparison of PT with the skewness 
estimated in the mock maps made with the weak lensing 
efficiency function. The errors are larger here than in the 
galaxy case because the depth (and volume) is larger and 



6 



Skewness of the projected density 




0(deg) 0(deg) 



Fig. 4. Angular skewness S3 (bottom) and kurtosis S4 (top) in mock maps with the APM power spectrum (left panels) and 
the CDM power spectrum (right panels) for the galaxy catalogue selection function. Results are symbols with errorbars and 
are compared with small angle approximation PT predictions (continuous line). The short-dashed line shows the Monte Carlo 
integration. The long-dashed line are results from smaller mock catalogues with the (APM) observational mask. 



therefore a given physical scale corresponds to a smaller 
angular scale. 

Sampling variance 

We find that it is crucial to use a large number of cata- 
logues in order to have a robust estimation at the largest 
angular scales. This is illustrated in Figure [| which shows 
moments of the distribution of relative errors in S3, i.e. 
AS3 = S3/S3 — I., from map to map around the mean 
value, S3. The rms error (dashed line) increases with an- 
gular scale as expected. The skewness of the error distri- 
1/3 

bution, ^(As 3 ) 3 ) (points in the figure), has large fluc- 
tuations, but seems to have a tendency to go negative on 
the largest scales in each catalogue. We show results for 
both the large APM simulations (open circles and long 
dashed lines) and both CDM and APM L = 400ft" x Mpc 
simulations. This result indicates that it is more likely to 
find smaller values of s 3 , when doing a smaller sampling, 
e.g. in a single map. 

The above arguments are also illustrated by compar- 
ing the results from the full sky mock catalogues with the 
mean values using 10 mock catalogues with the APM mask 
(e.g. GB97). The latter not only cover a smaller fraction of 
the sky (only about 10%) but are also subject to bound- 
ary effects. The results for these realistic maps are shown 
as long-dashed lines in Figures ||. The finite volume and 
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Fig. 5. Angular skewness S3 for weak lensing efficiency func- 
tion from simulations compared with small angle approxima- 
tion PT predictions (continuous line). The short-dashed line 
shows the Monte Carlo integration. 
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boundary effects seem important on scales 6 > 1 deg for 
CDM and about 9>2deg for APM. 
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Fig. 6. Relative errors in the the distribution of values of s3 
from different maps as a function of the angular scale. Dashed 
line corresponds to the rms value: ((As-j) 2 ^ ^ , whereas the 

points correspond to a measure of the skewness: ^(As3) 3 ^ 3 ' . 
Open circles and long dashed lines correspond to the maps from 
the largest box simulations. Filled triangles and short-dashed 
lines show resulsts from the small boxes. 



4-2. Time evolution 

Perturbation theory predicts time independent values of 
the skewness, in the limit of small variance. For large vari- 
ance, N-body simulations find that the 3D skewness, S3, 
increases with time due to non-linear effects (e.g. Baugh, 
Gaztanaga & Efstathiou 1994, Colombi et al. 1997). Non- 
linearities also erase the shape dependence in the hierar- 
chical structure, e.g. the bispectrum, (see Scoccimarro et 
al. 1997) reproducing a simple tree-level hierarchy, inde- 
pendent of the shape of the configurations. When this hap- 
pens there is a simple relation between the angular and 3D 
skewness: S3 ~ R3 S3 (see G94), so that the corresponding 
2D projected amplitudes are closer to the 3D PT results. 
Thus, on quasi-linear scales, the relation S3 ~ R3 S3 typi- 
cally underestimates the projection effects, while it should 
be more accurate on smaller scales, were PT results will 
overpredict the projection effect. 
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Fig. 7. Angular skewness (bottom) and kurtosis (top) in CDM 
mock maps from different output times: as = 1 (filled cir- 
cles) and as = 0.5 (open triangles). Continuous lines show the 
small angle approximation PT predictions. The short-dashed 
line shows the Monte Carlo integration. 

Figure [7] shows a comparison of S3 and S4 from mock 
maps drawn from two different outputs in the CDM galaxy 
model. Contrary to what one would expect, the later out- 
put seems closer to PT results. This is because of the pro- 
jection effects mentioned above, which compete with the 
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non- linear growth in S3. In fact, the slight disagreement at 
the smallest angular scales between the results of the full 
sky map (symbols with errorbars) and the one with the 
APM mask (long-dashed lines) shown in Figure |] is also 
due to the fact that we have used slightly different output 
times in each case. 



5. Conclusion 

We conclude from this analysis that the P.T. results for 
the skewness and kurtosis are accurate at degree scales 
or larger for the projected density. Comparisons between 
Monte-Carlo integrations and results obtained from the 
small-angle approximation are in good agreement. A de- 
parture starts to be noticeable only for smoothing scales 
above 5 degrees, but even there it is small compared with 
typical errors. 

The iV-body results allow us to estimate the small an- 
gular scales at which nonlinear effects start to affect the 
values of S3. We found these effects to be important below 
0.5 deg for the lensing case, and below 1 deg for the galaxy 
selection function, with a rapid growth of S3 at small scale. 
This however depends strongly on the initial power spec- 
trum. This effect is found to be more important for the 
measured APM power spectrum than for the CDM one. 

We have found two effects that might prove important 
in comparing perturbation theory with observations of an- 
gular clustering, and in particular for the APM. First, vol- 
ume and boundary effects are important on scales > 2 deg 
and tend to produce smaller values of S3 and S4 . We argue 
that this is not because the mean is biased, but because the 
distribution of errors seems to be negatively skewed (e.g. 
Figure ^) and it is therefore more probable to find smaller 
values. Second, the tree-level hierarchy model for projec- 
tions commonly used in the literature (e.g. by Groth & Pe- 
bles 1977, Fry & Peebles 1978, Szapudi, Szalay & Boschan 
1992, Szapudi et al. 1996), and in particular for the APM 
(G94) is not accurate on quasi-linear scales, as indicated 
in B95, because it underestimates the projection factors. 
At small scales, non-linear effects tend to mix these two 
cases (e.g. Figure]]]), while at larger scales this competes 
with large volume effects. The errors involved are of the 
order of 20% in S3, but more quantitative analysis of this 
point and implications for the APM observations will be 
presented elsewhere. 

To have a better insight into these nonlinear effects, 
and how it depends on the shape of the power spectrum, 
it could be fruitful to extend the recent results of Scocci- 
marro (1997) on the one-loop correction of the bispectrum 
to the projected third moment. 
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Appendix: The expression of the three-point cor- 
relation function in real space 

We want to express the three point correlation func- 
tion provided by Perturbation Theory (Peebles 1980, Fry 
1984), 



£ 3 (ri,r 2 ,r 3 ) 



d 3 ki d 3 ki 



(2tt) 3 (2tt) 3 



(25) 
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where 

Xl2 • Xi 3 

Xl 2 X 13 



u = 



(35) 



in real space only. The reason is that the integration in k 
space with real space top-hat window function is hardly 
possible because the integrals converge very slowly. 
Let us introduce the function <fi(r), 



f d 3 k P{k) 
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Then we can notice that 
d 3 k 
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VC(x) = -i 
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x dj{x) 
x dx 



k exp(ik • x) 



(and a similar property for <p) so that 
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This result generalizes the one obtained by Fry (1984) for 
power law spectra. 

The Monte-Carlo computations of the geometrical av- 
erages of £ 2 and £3 can then be done in real space, provided 
the functions £ 2 , £ 2 an d f' are known. Then the integra- 
tions can be reduced to a 6-dimensional integrals by direct 
integration over the azimuthal angles (see fig. 3). 



and one recognizes one term that intervenes in the expres- 
sion of the three-point function. 

To complete the calculation one can also notice that, 



Vij( x 12)<Py( x 13) 
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Finally remarking that 
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Fig. 8. Description of the variables used in the Monte Carlo 
Integration. The distance of xi to the origin, the distance be- 
tween xi and X2, the angle *I>i and the angle are all chosen 
randomly. The integration over <pi2 can be done explicitly since 
the distances of X2 to the origin and to xi remain then con- 
stant. 



one can express the three point correlation function in real 
space with the function £2(2), £ 2 ( x ) and ip'(x) only: 

£ 3 (ri,r 2 ,r 3 ) = y£ 2 (x 12 ) £2(0:13) + 

4 ip'{x 12 ) (fi'{x 13 ) 
7 X12 xi 3 



y (6(^12) + 2 



<p'(x 



i2, 



Xl2 



£2(0:13) + 2 



(fi'(xi 3 ) 



X13 



-u [£ 2 (x 12 ) <p'{x 13 ) + C 2 (xi 3 ) <p'(x 12 )] + 



6(^12) + 3 



<f'(x 12 ) 
X12 



£2(0:13) + 3 



<P'(X13) 
X13 



(34) 

This article was processed by the author using Springer- Verlag 
C y c ->14T E X A&A style file L-AA version 3. 



